Width difference in tfie B^ — B^ system from lattice HQET 
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We present recent results for the prediction of the B^ — B^ lifetime difference from lattice Heavy Quark Effective 
Theory simulations. In order to get a next-to-leading order result we have calculated the matching between QCD 
and HQET and the two loop anomalous dimension in the HQET for all the AB — 2 operators, in particular 



for the operators which enter in the width difference. 
(5.1 ±1.9 ±1.7) X 10"^ 



We obtain for the B^ — B^ lifetime difference. 



1. Introduction. 



'AFb, 



The width difference (^^) in the B° - 



B2 



system is expected to be the largest among bot- 
tom hadrons. The Standard Model prediction 



for 



AFe 



relies on an operator product expan- 



sion, where the short distance scale is provided 
by the b quark mass jl|. The theoretical esti- 
are in the range 5 -;- 15% and 



AFe 



mations for 

hadronic matrix elements of four quark operators 
are crucial inputs. We present a next-to-leading 
order (NLO) computation of the dominant matrix 
elements which contribute to -n-^. 



2. Standard Model formula for the width 
difTerence in the B'^ — B'^ system. 

The theoretical expression for the width differ- 
ence reads [|] 
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G{z)^B{mt) + Gs{z)^Bs{mh)X + Si/^^ , (1) 



where 



X = 



Ml 



{nib + rUsY 



(2) 



*Govierno Vasco fellowship. 
TTalk presented by J. Reyes. 



z = ml/ml, G{z) and Gs{z) are NLO Wilson 
coefficients calculated in [g|, ^i/mj are 0{l/mf,) 
corrections which depend on matrix elements of 
higher dimension operators (the explicit expres- 
sion for 5i/m^ can be found in [Q). 
B and Bs are bag parameters defined as vacuum 
insertion deviations, 



{Bs\OLi^i)\Bs) EE -flMlBifi), 



= -^fsM^X Bsifi) 



{Bs\Osi^i)\B. 

where 

Ol = 6''7"(l-75)g*fe^'7M(l-75)'?^' 



(3) 
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i, j are color indices. 

One of the most important sources of uncer- 
tainty in (|l|) is the Bg decay constant, fs^- In 
order to reduce it we follow ref. 101 and write 



AFe 



ATb: 



in the following way. 
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iG{z) + Gs{z)'R{mb) + ~5i/r, 



(5) 



Notice that all the non-perturbative contribution, 
at leading order in l/nib, comes from the TZ pa- 
rameter. 



TZimb) 
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X, (6) 
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= , is accurately com- 



The parameter ^ = 



puted from lattice simulations Q, a;^ = ml/m^ 
and, rj^Omi,) and So{xt) are well known theoreti- 
cally [glJG). Finally AniBa and tb^ are accurately 
measured experimentally. 

3. NLO matching from continuum QCD to 
lattice HQET. 

Our aim is to compute the two B-parameters 
B and Bs from lattice simulations. The b quark, 
however, is too heavy to be simulated dynami- 
cally in accessible simulations since the typical in- 
verse lattice spacings, are in the range 2-^3 Gev. 
There are two ways to overcome this problem: 
performing simulations around the charm mass 
and extrapolate the results to the physical b 
quark mass or using an effective field theory 
(HQET, NRQCD) in which the heavy degrees of 
freedom have been integrated out. In our work 
we follow the second method using HQET. Our 
task is to extract matrix elements of operators in 
QCD from lattice HQET simulations. To obtain 
a NLO result the following four steps, described 
in detail in H, should be computed: 

1. The continuum QCD-HQET matching. 

The QCD operators are expressed as linear 
combinations of HQET ones in the contin- 
uum at a given high scale, say, fj, — rrib. 

2. Running down to /i = a ^ in the HQET. 

The HQET operators obtained in step 1 at 
the high scale /i = mt are evolved down to a 
lower scale ^ — a^^, appropriate for lattice 
simulations, using the HQET NLO renor- 
malization group equations, which implies 
the computation of the two loop anomalous 
dimension of the relevant operators. 

3. Continuum lattice HQET matching. 

Having obtained the continuum HQET 



operators at the scale fx 



they are 



expressed as a linear combination of lattice 
HQET operators at this scale. 

4. Lattice computation of the matrix 

elements. 

The matrix elements of the lattice HQET 



operators in Step 3, are measured by Monte 
Carlo numerical simulations on the lattice. 

We stress that the first three steps depend on 
the renormalization scheme used in the HQET. 
When the two loop anomalous dimensions are 
included in the renormalization group evolution 
(RGE) this scheme dependence must cancel up 
to 0{as)- So far, the complete chain of match- 
ing equations have been calculated for Ol only. 
However, the two last steps are known for all op- 
erators y,^,^. We have extended the calculation 
to all the AB — 2 operators computing the con- 
tinuum QCD - HQET matching and the two loop 
anomalous dimensions |l^,fl| . Wc have explicitly 
verified the cancellation of the dependence of our 
expressions on the intermediate HQET scheme 
]10| . This is a strong check of our calculation 
of the matching. 

4. B-parameters from numerical simula- 
tions on the lattice. 

The lattice B-parameters are extracted from 
the large time behavior of the ratio between three- 
and two-point correlation functions (see for exam- 
pie i): 
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where 



Co.(ii,i2)=5I<0|v4o(a?i,ii)0.(0,0)4(x2,t2)|0) 
C{t) = Y.(^\Ao{x,t)Al{0,0)\0} 
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and A^ is the HQET axial current. 

The results presented here are obtained from a 
quenched simulation performed by APE collabo- 
ration in a 2^ x 40 lattice with 600 gauge con- 
figuration at /3 = 6.0 with a^^ ~ 2 Gev. The 
details of the simulation can be found in |j] and 
references therein. 

An important remark we want to make is 
that, as firstly pointed out in Q, all lattice B- 
parameters obtained in the static limit are very 
close to their vacuum insertion values. This re- 
sult has been obtained in static simulations by 



APE coll. H, UKQCD coll. |j, Christensen 
et al. |l3| and also by JLQCD coll. in NRQCD 
when their results are extrapolated to the infi- 
nite mass limit |1J]. Notice that although these 
simulations are performed with different actions 
(static-clover, static- Wilson, NRQCD-clover) at 
different lattice spacings and by different groups, 
the conclusion is the same in all cases: the lattice 
B-parameters in the static limit are compatible 
with the vacuum insertion value within the sta- 
tistical errors. This surprising property of the 
HQET deserves further investigation which is in 
progress. 

5. l/nib dependence of the B-parameters. 

Having performed the calculation described in 
the previous section one obtains the B-parameters 
in QCD to leading order in 1/mf,. By using 
vacuum insertion approximation (VIA) for the 
subleading operators |l^] and the fact that the 
bare lattice B-parameters are very close to the 
VIA value, we have shown iQ that all AB = 2 
B-parameters, defined as vacuum insertion de- 
viations, have small 0[A/mb) corrections (A = 
Mb, -mi,), 

B,{mb) ^B,{mb) + 0(0.3 — )+0{^), (7) 

^ ^ ■> ^ 

~0.05 

where Bi{mb) and Bi{mb) are the full and the 
static QCD B-parameters respectively. The fac- 
tor 0.3 is an estimation of the deviation from VIA 
in the calculation of the subleading operators and 
0{as) corrections. 

6. Results. 

Before presenting our results, let us stress that 
rris and ttij, parameters in previous sections are 
the corresponding pole masses. The latter coin- 
cides with the expansion parameter of the HQET 
because we have set the residual mass to zero. We 
calculate the b quark pole mass from the running 
MS mass, which can be accurately determined. 
Since our computation is performed at NLO, we 
use the perturbative relation between the pole 
and the running mass at the same order. From 



Table 1 

Values of the B-parameters in the NDR-MS 
scheme (see text). 

B-parameters 



B{mb) 
Bs{mb) 
B^'^imb) 



0.83(5)(4)(5) 
0.96(8)(5)(5) 
0.94(5)(4)(5) 



g|^(mb) 1.03(3)(5)(5) 



the world average running mass [Q, mb{mb) = 
4.23 ± 0.07 Gev one obtains mb = 4.6 ± 0.1 Gev. 
Notice that the contribution of tUs is very small 
because it always appears divided by mb- 

In table |^ we present the values for all AB — 2 
B-parameters. We have included, for the sake of 
completeness, the other AB = 2 operators in the 
HQET, 



Bs\0'^''if,)\B, 



^n-^ 



^ S^''(A^),(8) 



(^4op4^ 



-1 _ oLR 



Ifk^k 



srip^). (9) 



and 



qLB. ^ l^^^^^l_^^)q^y^^^l^^^)qJ^ (10) 
^LR _ ui ^-1 „, \„iu n I _, \ „.i Qi') 



or = 6^(l-75)<7*6^(l+75)<7^ 



The first error in table 1 comes from lattice sim- 
ulations and includes statistical and systematic 
errors. The second one is an estimate of the error 
due to the uncertainties in the values of the lattice 
coupling constant and to higher order contribu- 
tions to the matching. The third is an estimate of 
l/wf, corrections to the static result (see eq.(7)). 
The analysis of the perturbative matching is ex- 
plained in detail in our previous work |^ where 
we gave B{mb) = 0.81 ± 0.05 ± 0.04. The tiny 
difference is due to the fact that there we used, in 
the perturbative evolution, a number of fiavours 
Hf — A instead of nj = as in the present paper. 
The B-parameters in table || are QCD scheme de- 
pendent, they all are in the NDR-MS scheme. In 
order to make contact with the NLO computation 
of Q, we have subtracted the evanescents opera- 
tors in the renormalization of O and Os as in p|. 
For B^^ and Bg^ we use the prescription of |17|. 



Table 2 

Values of B, Bs and TZ from different groups. 



Groiip B^n'Hi) 



Bsjrnb) 



Tlimt) 



[]|4]J 0.85(3)(11) 0.95(2)(12) -0.91(5)(17) 



3| 0.91(3)+°:°° 1.02(2)i°:°l 



-0.93(3)i°:°° 



This work 0.83(5)(6) 0.96(8)(7) 



-0.95(7)(9) 



For the parameter TZ (see eq.(6)), which con- 
tains all the non-perturbative contribution to the 
width difference, at leading order in 1/mi,, we ob- 
tain 



n{mb) = -0.95(7)(9) 



(12) 



where as usual the first error comes from lattice 
simulations, and the second is systematic due to 
the uncertainty in the perturbative matching and 
to the contribution of higher orders in l/rrif,. 

6.1. Comparison with other recent results 

In table y we present our result for the B- 
parameters relevant for „ "°° compared with 

other recent determinations. In Q the Bs is de- 
fined as in eq. (|j) but in terms of the running mass 
instead of the pole mass. Therefore, we have mul- 
tiplied their value by {mb/mi,{rn(,))'^ in order to 
compare with our result. On the other hand, in 



the definition used in |14| does not appear the 
factor X (see eq. @)), therefore we have divided 
their result by X. Notice that in spite of using dif- 
ferent methods to obtain the B-parameters, there 
is a good agreement between the three computa- 
tions, in particular in the value of the ratio TZ. 
From equation (0) we get our prediction, 

Ar 

-— ^ = (5.1 ±1.9 ±1.7) X 10^2. (13) 

The first error is systematic obtained from the 
spread of values of all input parameters in eq.(||). 
The second one comes from the uncertainty in 
the value of (5i/mf Since in the estimate of ^i/mj 
the operator matrix elements were computed us- 
ing VIA and the radiative corrections were not 
included, we assume an error of 30% Q. As can 
be seen, this parameter is still affected by a large 
uncertainty, so that a precise determination of the 
width difference requires the computation of the 



subleading matrix elements using lattice QCD. 
This simulation is missing to date. 

Our result is to be compared with the present 
experimental status ||l8|| 
exp. 

(14) 



AF, 
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l7tTo) X 10- 



As can be seen the central values are rather dif- 
ferent but still compatible within the large errors. 
More work is needed on the theoretical and ex- 
perimental sides to reduce the uncertainty in the 
width diference in the Bs system. 
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